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GREEN’S FUNCTIONS FOR 
STIELTJES BOUNDARY PROBLEMS 

M. ROSENKRANZ AND N. SERWA 


Abstract. Stieltjes boundary problems generalize the customary 
class of well-posed two-point boundary value problems in three 
independent directions, regarding the specification of the boundary 
conditions: (1) They allow more than two evaluation points. (2) 
They allow derivatives of arbitrary order. (3) Global terms in 
the form of definite integrals are allowed. Assuming the Stieltjes 
boundary problem is regular (a unique solution exists for every 
forcing function), there are symbolic methods for computing the 
associated Green’s operator. 

In the classical case of well-posed two-point boundary value 
problems, it is known how to transform the Green’s operator into 
the so-called Green’s function, the representation usually preferred 
by physicists and engineers. In this paper we extend this trans¬ 
formation to the whole class of Stieltjes boundary problems. It 
turns out that the extension (1) leads to more case distinction, (2) 
implies ill-posed problems and hence distributional terms, (3) has 
apparently no effect on the structure of the Green’s function. 


1. Introduction 

Boundary problems for linear ordinary differential equations (LODEs) 
or linear partial differential equations (LPDEs) are certainly among the 
most important model types in the engineering sciences. Interestingly, 
their systematic treatment in Symbolic Computation started rather 
recently m- For handling the central problems of solving and fac¬ 
toring boundary problems, a differential algebra setting for LODEs is 
employed in [201 HH] and for LPDEs in [211 CH]- An overarching ab¬ 
stract framework based only on Linear Algebra is developed in [15] . 
For the classical treatment of boundary problems in Analysis, we refer 

to [niiniioiiioii. 

In this paper we restrict ourselves to LODEs, where the “industrial 
standard” for solving boundary problems is their so-called Green’s func¬ 
tion. This is in stark contrast to the operator-based methodology used 
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in the above references. In fact, given a fundamental system, the al¬ 
gorithm of [m EH] computes the solution of a boundary problem in 
the form of its Green’s operator. In the classical setting of well-posed 
two-point boundary value problems (see Section [2]), this algorithm ad¬ 
mits an optional extra step for extracting the corresponding Green’s 
function. Our goal here is to extend this postprocessing step to the con¬ 
siderably larger class of Stieltjes boundary problems (see Section [2]). 

One way to understand the relationship between Green’s operators 
and functions is to view the latter as a certain canonical form. For mak¬ 
ing this precise we equip the ring of integro-differential operators with 
a slightly different set of reduction rules favoring multiply initialized 
integrals, leading to the ring of equitable integro-differential operators 
(Section [3]). 

A simple example will make this clear—in fact the simplest of all 
honest boundary problems [13 §3.2]. Given a forcing function / G 
G°°[0,1], we want to hnd u G C'°°[0,1] such that 

u" = /, 

m(0) =’m(1) = 0. 

The Green’s operator G: C'°°[0,1] —)■ (7°°[0,1] of this problem is de- 
hned by Gf = u. Using the standard reduction system of pni , we 
would distinguish one integral like jf := and then obtain 

the Green’s operator in the canonical form 

(1) G = - jx + x[l\^x-x\_l\j, 

where [aj denotes the evaluation functional / i—)■ /(a) for any a G R, 
in analogy to the multiplier notation of nn. For extracting the Green’s 
function, however, it is more useful to use the alternative canonical form 

( 2 ) G = xi^x-xJ^x-^^x + xj^ 

where Lf ■■= j : J{f)df now denotes the integral initialized at the 
point a G {0,1}. In fact, this is the form given in [T7], and we shall 
see in Section |3] that the setting of biintegro-differential operators used 
there is essentially a special case of the equitable operator ring em¬ 
ployed in this paper. The point of the canonical form ([2]) is that it 
allows us to apply the dehning relation Gf{x) = df of 

the Green’s function directly to obtain the latter as 


{x - 1 )^ 
a;(^ - 1) 


for 0 < .^ < a; < 1, 
for 0 < X < .^ < 1. 


( 3 ) 
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Heuristically speaking, one moves the terms to the upper and the 
terms to the lower branch, at the same time translating the “a:” after 
the integrals into Note incidentally that g{x,^) = g{^,x) in the 
above Green’s function ([3]). As is well known in Analysis [5l §7] [2ll 
§5], this is a consequence of the self-adjoint nature of this boundary 
problem—a topic that we would wish to investigate in the future for 
the more general class of Stieltjes boundary Problems (Section [6]). 

We will elaborate on the above principles to generalize it in three “or¬ 
thogonal directions”: (1) We allow more than two evaluation points, 
leading to an increased number of case branches. (2) Using deriva¬ 
tives of arbitrary order in the boundary conditions leads to distribu¬ 
tional terms. (3) Boundary conditions with integral terms (so-called 
“non-local problems”) are also included; they do not lead to further 
complications. 


2. Stieltjes Boundary Problems 

For giving a precise dehnition of the class of admissible boundary 
problems, we follow the setting of [20] . Hence let {iF, d, J) be a fixed 
ordinary integro-differential iP-algebra (here ordinary means kerd = 
K). Later we shall specialize this to = G°°(1R), the real- or complex¬ 
valued smooth function. This is theoretically convenient but of course 
needs to be replaced by a suitable constructive subalgebra for actual 
computations. 

The ring of integro-differential operators over T, introduced in [201 
§3], will be denoted here by J] to emphasize its dependence of 

the chosen set of characters <F, and also to mark the contrast to the 
equitable operator ring J^[d, to be introduced in Section [3l where 
the integral operators are parametrized by <h. In the case of = 
G°°(R), these characters will be evaluations at given points of IR so 
that we may take $ C ]R. 

We recall the standard decomposition 

(4) J-,[c),Jl=J^[81+JF|J] + (<l.), 

where J^[d] denotes the subalgebra of differential operators (the K- 
subalgebra of J] generated by and 9), J^[j] the nonunital subal¬ 
gebra of integral operators (the nonunital A'-subalgebra of J] gen¬ 
erated by T and J), and (d*) the two-sided ideal of J] generated 
by the characters in $. The corresponding right ideal |<I>) = <!>■ J] 
is known as the ideal of Stieltjes conditions^ and one may check that (<I>) 
is in fact the left J^-module generated by the Stieltjes conditions. 

From the viewpoint of applications, Stieltes conditions (3 G ($) are 
easier to comprehend in terms of their J]-normal form: They can 
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be described uniquely as sums 


(5) /3 = + 


with only finitely many G K and nonzero. The double sum 

in ([5]) is called the local part of /9, the subsequent sum its global part. 
In the important C°°(1R) case with distinguished integral J = this 
yields 

(p^i p 

for certain E IR and E C'°°(1R). 

An n-th order Stieltjes boundary problem is a pair (T, B) with a monic 
differential operator T G J^[d\ of order n and a bonndary space B < T* 
given as linear span B = [(3\,, /5„] of n linearly independent Stieltjes 
conditions. In traditional representation, snch a bonndary problem is 
displayed as 


Tu = f 
jSiU = ■ ■ ■ = (3nU = 0 


with the nnderstanding that u G -F is desired for any prescribed forc¬ 
ing fnnction f E J^. For the (usnal) Green’s operator to be well- 
defined, we need the bonndary problem ([6]) to be regular in the sense 
that kerT -j- B^ = where B-^ = {u E \ f3{u) = 0 for all (3 E B} is 
the corresponding space of admissible fnnctions. Regnlarity is eqniv- 
alent to the reqnirement that ([6]) has a nnique solution u G -F for 
every given f E 3F. An algorithmic method for testing regnlarity starts 
from a fnndamental system ui, ... ,Un E 3F for T, meaning a iF-basis 
of ker T. Then ([2]) is regular iff the evaluation matrix 


(7) 



(3l {Un) 


fdn ig^n) 


^ j^nxn 


is regnlar; see (15) of [20]. Given a fnndamental system ui,...,Un 
for T, the solntion algorithm of [20] compntes the Green’s operator of 
any regular Stieltjes boundary problem as an integro-differential oper¬ 
ator G E [9, J]. 

Within the class of Stieltjes bonndary problems, we make the fol¬ 
lowing distinctions in order to characterize the classical scenario as a 
certain special case. 
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Definition 1. A Stieltjes boundary problem (T, B) of order n with B = 

, /9„] is called well-posed if the (di can be chosen with all deriva¬ 
tives having order below n; otherwise it is called ill-posed. Furthermore, 
we call (T, B) an m-point boundary problem if the maximal number of 
evaluation points occurring in any iF-basis (/?,) of B is m, and we 
call (T, B) local if the fdi can be chosen without global parts. 

Let us digress a bit on the notion of ill-posed boundary problems. 
Following Hadamard, a problem is generally called well-posed 0 p. 86] 
if it is regular (meaning its solution u exists and is unique for all given 
data /) as well as stable (meaning u depends continuously on /). Oth¬ 
erwise, one speaks of an ill-posed problem. In the case of boundary 
problems ([6]), we search for u G C'°°(]R), and the data is given by the 
forcing function / G 0°°(]R). Stability—and hence well-posedness— 
depends on the topology chosen for the function space C'°°(]R). Using 
the norm as in many application problems, the distinction between 
well- and ill-posed boundary problems coincides with the one given 
above. 

Since local boundary problems involve only evaluations of the un¬ 
known function (rather than definite integrals), we also call them “bound¬ 
ary value problems”. We can now characterize the classical case, de¬ 
scribed for example in 0 §7], by the following three-fold restriction: 
They are the well-posed two- point boundary value problems. (Some¬ 
times one meets the further restriction to self-adjoint boundary prob¬ 
lems.) 

The classical case (in the above sense) is clearly the most frequent 
case in the applications (but this could also be due to a selection bias: 
having a well-equipped toolbox for classical problems might tempt en¬ 
gineers to restrict their attention to classical problems). Nevertheless, 
multi-point boundary value problems are also important for some ap¬ 
plications lU Ex. 1.6], [23], [13], [H], |3]. Boundary problems with 
nonlocal conditions are more seldom, they are usually studied for non¬ 
linear equations H, mi; the linear case serves as the initial approxi¬ 
mation. Finally, the case of ill-posed boundary problems is—for obvi¬ 
ous reasons—mostly avoided when engineering problems are modelled. 
However, there are cases where their treatment is inevitable, typically 
in the context of inverse problems 0 . Since the numerical treatment of 
such problems is very delicate, it is of paramount importance to have 
exact symbolic algorithms wherever this is possible. 

We will lift all three of these restrictions for the algorithm of ex¬ 
tracting Green’s functions, which will be given below (Section Hj). As 
indicated in the Introduction, the crucial tool for this purpose—even in 
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the classical case—is the ring of equitable integro-differential operators 
with its alternative canonical forms. 

3. Equitable Operators 

The passage from the standard integro-differential operator ring J] 
to its equitable variant J^[d, is based on the fundamental theorem 
of calculus /^(O= /(^) “ f{y^) for function / G 0°°(]R) 
and initialization point (p G R. Likewise, if {T, d, J) is an arbitary 
integro-differential iL-algebra and ip a character (multiplicative linear 
functional), one can use the dehnition := (id—(p)J to obtain the 
corresponding relation 5 = id —p. In some contexts (especially in 
the presence of several integral operators like and on bivariate 
functions), it may be useful to write the integral as If is 
another character, one observes the relation — Sp ~ ^rid it is 
natural to write for both expressions. 

Note that (J^, d, J^) is again an ordinary integro-differential iL-algebra, 
and the preference of J over can appear arbitrary in certain settings. 
Accordingly, one may build the ring of integro-differential operators by 
adjoining all while the characters p themselves are now redundant 
due to the above fundamental relation. The precise formulation of the 
resulting ring J^[d, as a quotient is described in |T6l §5.1]. For our 
present purposes, we shall only list its relations (see Table [T] where • 
denotes the natural action of the operators), which are an easy con¬ 
sequence of the relations of the standard integro-differential operator 
ring Jj. 


fg f -g 

df ^ d-f + fd 

5 j : ^ 1 

rX n rX / rX 

Jip ^ JTp ^ U ip 

/) j; - j : ( j : • /) 


rjd ^ f-rjd.f)-p.f+{p-f)f:;d 


Table 1. Equitable Integro-Differential Relations 


Similar to the standard decomposition (|1]), we have also the equi¬ 
table decomposition = J^[d] + where is 

the nonunital subalgebra of equitable integral operators Yl'i=o 
and the J^-submodule consisting of 'f2i=o /*leads to 

the obvious normal forms in J^[d, J^j. 
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The so-called translation isomorphism l: J] —)■ leaves 

f E J- and d invariant while using the above fundamental relation 
via L{ip) = id — d and l ^(J^) = (id— 99)/. Note that this holds 
also for the character € := id—associated to the distinguished 
integral J = Xg underlying [5, J]. 

Specializing to = C°°(1R) and <h = {0,1}, we can deal with the ex¬ 
ample in the Introduction, where we have the normal form G G iF<s>[d, J] 
in ([1]) along with its equitable variant l{G) G J^] in ([2]). In 
such two-point cases with characters $ = the equitable op¬ 

erator ring J^[d, J,j,] is essentially the same as the ring of biintegro- 
differential operators. More precisely, we obtain a biintegro-differential 
algebra (J^, d, J*, J^) with integral j* '■= and cointegral := — 

in the sense of [161 Def. 3.23]. Note that J* and are adjoint with 

respect to the inner product dehned by {f\g) := (J* +J*)(/ 5 ') = Stfd- 
Incidentally, the notion of biintegro-differential algebra coincides with 
the (badly named) notion of “analytic algebra” introduced in [13 
Def. 2] and replicated in [201 Ex. 5]. Clearly, the operator ring re¬ 
sulting from J^[d] by adjoing J* and is the same as J^[d, J^j, modulo 
the sign change in the cointegral. 

Note that also G J] and ip G J^[d, are legitimate expres¬ 

sions via the above translation isomorphism. They are not in canonical 
form but we may think of them as a kind of abbreviation for the corre¬ 
sponding canonical expression. 

In fact, the extraction of Green’s functions is based on the fol¬ 
lowing slight variation of the equitable integro-differential operator 
ring -^[5, J^j. Writing any element U G -^$[5, J] in the form U = 
T + K + B with T G -T[9], K G -F[J] and B G (<h) according to (jl]), 
we let T G J,j,] and K G C invariant while translat¬ 

ing B G (<h) as follows. Since (<h) is the left J^-module generated by 
Stieltjes conditions (jSD, we may split B = A-I -7 into a left J^-linear com¬ 
bination A of local Stieltjes conditions and a left J^-linear combination 7 
of global Stieltjes conditions. It turns out to be expedient to keep A 
in this form, without eliminating the characters via p = id—but 
to translate 7 via ='■ This is what we mean when 

referring in the sequel to the equitable form of an integro-differential 
operator U. 


4. Extracting Green’s Functions 

We now turn to the central task of this paper, the extraction of the 
Green’s function g{x,^) corresponding to the Green’s operator G G 
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J^[d, J#] computed by the algorithm of [ 20 ] and converted to equitable 
form as described in Section |3l Hence we specialize now to = C'°°(]R). 
Note that we may think of g{x,^) as a kind of coordinate representation 
of the induced operator action G: iF ^ iF] in quantum mechanics this 
would correspond to the “position basis” (as opposed to the “momen¬ 
tum basis” in the Pontryagin dual reached via the Fourier transform). 
Hence we will use the notation g{x, = Gx^, thinking of the x, ^ rather 
like continuous indices similar to the discrete indices i,j in the matrix 
elements Aij of some A G 

In fact, we will use this notation Gx^ for any equitable integro- 
differential operator G G J^]. Its result will in general contain 
Dirac distributions [22l §2] and their derivatives but nothing beyond 
that. Since all boundary problems considered in this paper have only 
hnitely many evaluation points a G <F C R, one may choose an inter¬ 
val J C ]R containing all the a. Hence the C(J^)-module Q C 
generated by the Dirac distributions 6a and their derivatives will be 
sufficient to capture all Green’s “functions” Gx^ G Q. Here and in 
the sequel we shall follow the common engineering (and also applied 
maths) practice of referring to distributions like 6a as functions. In the 
same vein, we shall also write 6{^ — a) in place of 6a, in view of the 
dehning property Jj ^(^ - a) /(O d^ = f{a). 

The transformation from Green’s operators to Green’s functions 

G, G Gx^ 

is clearly an R-linear map, hence it will be sufficient to dehne it on the 
canonical ^-hasis of 6F^[d, J] = 6F[d, J^]. Following the strategy of the 
example in the Introduction, the easiest part is -F[J] C 6F<^[d, J], which 
is handled by setting 

U^9)xi = f{x) g{0 [0 < ^ < x] - /(x) g{^) [x < ^ < 0], 

where we use the Iverson bracket notation [P] signifying 1 if the prop¬ 
erty P is true and zero otherwise. Note that at most one of the two 
summands above is nonzero for hxed (x, ^). Since ($) C J] is a 

left J^-module over |<h), we settle this part via 

(/ L“JJs)iS = sgn(a) /7)sK) [0 < ^ < a]. 

Finally, on P[d] we dehne 

and the dehnition is complete in view of (jl]). Moreover, it is easy to 
check that the assignment G ha Gx^ is correct in the sense that Gf = 
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jjGx^ fiO The isomorphism l of Section [3] may now be employed 
to obtain the required transformation J^[d, —)■ Q. In fact, the above 
case fjgE J^[j] generalizes immediately to 

ifL9)x^ = f{x)g{0 [a < ^ < a;] - f{x)g{0 < ^ < a], 

which will turn out to be the essential clause for extracting Green’s 
functions of (well-posed) multi-point boundary problems. For seeing 
this, we need a more detailed description of the underlying Green’s 
operators. 

We turn hrst to the easy case of a one-point boundary problem, more 
appropriately known under the name of initial value problems (T, [€,..., €d"' 
for T G J^[d\ of order n. The corresponding Green’s operator is 
called the fundamental right inverse and can be computed eas¬ 
ily via the well-known “variation of constants” formula [T6l Thm. 6.4]: 

If Ml,..., is a fundamental system for T with Wronskian matrix IT, 
the fundamental right inverse is given by 




u 


J=1 


d' 


Here d = det(lT) and di = det(lT), where Wi denotes the matrix 
resulting from W when replacing the Gth column by the n-th unit 
vector of K^. 

What we shall need in the sequel is how reacts to left multipli¬ 
cation by J^[d]. 


Lemma 2. LetT G iF[d] he any monic differential operator of order n, 
and choose a fundamental system Ui,... ,Un for T with Wronskian ma¬ 
trix W. Then we have 


( 8 ) 




d 


with 


Pk 


1 


i=i 


(fc-i) 


dj G -F, 


where d and di are as above. 


Note that pi = • • • = Pn-i = 0 by the dehnition of the dj] hence the 
second sum in is only present for k > n, and we may equivalently 

write its range as j = n,..., k. Furthermore, we have pn = 1 from the 
definition of d. For k > n, however, the pk are functions of W, so in 
general they do not commute with the in the second summand 
of dHl). 
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Proof. We use induction on k. In the base case k = 0, this is the usual 
variation-of-constants formula as given in 0 p. 74]; see [201 Prop. 22] 
and |T 6 l Thm. 6.4] for its operator formulation. Now assume ([ 8 ]) for 
hxed A; > 0; we show it for k + 1. By the induction hypothesis we 
obtain 


Qk + lj.O ^ 

i=i 


^ •’ d 


+ 


d 


E 

i=i 


uf\lj 


+E9' 

i=i 


,k-j+l 


Pji 


which is just ([ 8 ]) for fc + 1 since the middle sum is pk+i and can be 
absorbed into the third. □ 


Lemma 3. The Green’s operator of any regular Stieltjes boundary prob¬ 
lem is contained in + C-, where C denotes the left iF-module gen¬ 
erated by the local Stieltjes conditions. 

Proof. Assume (T, B) is any regular Stieltjes boundary problem of or¬ 
der n with Green’s operator G, and let P be the projector onto kerT 
along B-^. By the proof of [201 Thm. 26] we have G = (1 — P)T^, 
and we know that P is an J^-linear combination of Stieltjes conditions 
by [201 (16)] in that same proof. From ([ 8 ]) it is clear that G -F[J^], 
so it suffices to show PT^ G -PlJ^] + C- Each summand of P is either 
of the form f[a\ or f[a\jg = fj^g - fj^g e P[J^]. In the lat¬ 
ter case we obtain an expression in P”!/^] since P[J^] is a (nonunital) 
subalgebra of P[cl, J^]. It remains to prove f[a\ G P[J^] + C. 

From dH]) we see that 

/bJ ^ fuf\a) /bJ d^~^Pr 

i=i i=i 

The first sum is clearly contained in P[J^], while the second is in C 
because pj G F[d] may be rewritten in canonical form as a sum of 
terms Pid’' so that [aj pj is a sum of local conditions a{gi) [aj 9* 
and hence itself local. □ 

We are now ready to state the main structure theorem for Green’s 
functions of regular Stieltjes boundary problems. 

Theorem 4. The Green’s function of any regular Stieltjes boundary 
problem with m evaluations Oi,..., has the form g{x,^) = g{x,^) -\- 
g{x,^), where the functional part g G G(P) is defined by the 2{m — 
1) case branches 

^ G bi, Oj+i] {fi<i <m),x <i] 
i G bi, Oj+i] (0 < i < m), ^ < X, 
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while the distributional part g{x, ^) is an J^-linear combination of the 8{f— 
Q!j) and their derivatives. 

Proof. If G is the Green’s operator of the given Stieltjes boundary 
problem, Lemma [3] says that G = G + G with G G and G E C. 

We will show that g{x,f) = G^^ and g{x,f) = G^^ are as described in 
the theorem. Starting with the former, we may write 

r 

i=l 

where Oj = aj is possible for i ^ j. Using the transformation 8F[d, J^] -E 
Q, we obtain g{x,f) as 

r 

X] (yft{x)gi{i) [ai < ^][^ < x] - fi{x) gi{f) < ai][x < 

i=l 

r 

= 9j{0) K-i < ^ 

2=1 aj<ai 
r 

< Oij+i][x < e], 

2=1 aj>ai 

where the two inner sums are restricted by j > 0 and j < n. Col¬ 
lecting terms, this is a sum of 2 (m — 1) characteristic functions over 
disjoint domains in R^, hence one may also write g{x,f) in terms of a 
corresponding case distinction with 2 (m — 1 ) branches. 

The distributional part g{x,f) is even easier. Writing G as an iF- 
linear combination of local conditions we obain g{x,f) via 

= (E = E(-l)‘- a). 

a,I ^ a,i 

which is clearly of the stated form. □ 

The above theorem is constructive, and we plan to implement the 
underlying algorithm on top of the Maple package IntDiffOp [T2] . 

Remark 5. If the distinguished character € = [Oj is not used in 
the boundary conditions, a straightforward translation of the Green’s 
operator G may introduce two spurious extra case branches in the 
Green’s function G^^ since occurs in the formula for G. For avoiding 
this, one has to use a different version of that replaces [Oj by any 
one of the characters [ajj used in the boundary conditions. 
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5. Examples 

Our first example (in addition to the minimal one from the Introduc¬ 
tion) is a four-point boundary value problem taken from [2], (2.1,2.2)], 
where we have specialized the parameters and rescaled the interval 
to J = [0,1] for the sake of simplicity. Hence we are dealing with the 
boundary problem 


-u” = /, 

n(0) + m(1/3) = m(1) + m(2/3) = 0, 


where we may assume u, f & C'°°[—2,2]. Note that this is a well- 
posed boundary problem, so the Green’s function will not have a dis¬ 
tributional part. Computing the Green’s operator with the IntDiff Op 
package yields after some rearrangements the result 

G = - P + (-5/24 + i/4) L1/3J J 

+ (5/8- 3a:/4)[l/3Jp + (1/8 - 3i/4)[lJJi; 

+ (1/12 - a:/2)[2/3J J + (-1/8 + 3i/4)[2/3jP 

Transforming G to equitable form is simple, via [crj J Jo-j: ;. We 
can then determine the corresponding Green’s function g{x, = g{x,^) 
with 6 cases, and its terms may be computed according to Theorem 01 
The result for g{x,^) is summarized in the table below. 


Case 

Term 

0 < C < 1/3T < X 

H 

1 

00 

0 < C < 1/3,2; < C 

(3/4)x^ + (3/8)C - X 

1/3 < C < 2/3,^ < X 

(3/2)x^ - (5/4)C - (l/4)x + 5/24 

1/3 < C < 2/3, X < C 

(3/2)x^ - (1/4)C - (5/4)x + 5/24 

2/3<C< l,C<a; 

(3/4)x^-(9/8)C + (l/4)x + l/8 

2/3 < C < 1,2; < C 

(3/4)xe - (1/8)C - (3/4)x + 1/8 


Our second example is, as it were, totally unclassical: It is ill-posed, 
has nonlocal conditions and contains three evaluation points —1,0,1. 
In our standard notation, we write this boundary problem as 

u" -u = /, 

- foU(^)^d^ = 0 , 

n'(-l)-n"(l) + /_\n(Ode = 0, 


where we assume now u, f G 1,1]. Using the method of [20] . 

it is straightforward to compute the Green’s operator G. In fact, the 
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IntDiffOp package yields the result 
aG = a/2(e^Je-^-e-^Je^) 

+ 2(-e"+3 + + 6-"+^ _ e-+i)([-lja + [Ij Jx) 

+ (e - l)(-e^+2 - 2e^+^ + e-^+^)(L-lJ J + UJ J) 

+ (3e"+2 _ e"+i - 3e-"+' + 36“") [Ij Je" 

+ (2e^+2 - 3e^+^)(e-H-lJ Je-^ + eL-lJ Je^) 

+ (-e^+3 - e^+2 + 2e^+^ + 6-^+^ _ g-^+i^ 

using the abbreviation a := 2(2e — 3)(e — 1) while collecting and fac¬ 
toring some terms for enhanced readability. After transforming this to 
equitable form (which is again straightforward), we can apply Theo¬ 
rem 0] to extract the Green’s function g{x,^) = g{x,^) + g{x,^) with 
the distributional part 

a g{x, 0 = (-e^+3 _ gX +2 ^ 2e^+i + 6-^+^ _ ^-x+i^ 

+ 2 (-e^+3 -F e^+2 - -F 6-^+^ - 6“^+^) 6 '{^ - 1) 

coming from the (...) [Ij and (...) |_1J9 terms, and with the functional 
part dehned by the case distinction for ag{x,^) as given in the table 
below. 


Case 

Term 

-1 < ? < o,c < * 

3g^+^+4; 3g^-« 2e“’+"-« - 2e’^+"’+« 

+e’*+"’ + 6"’"+^ + - 6“’’’+^ - 2e"’+^ 

-1 < ? < 0,® < ^ 

-2e"’+^ + 2e-“’+’'+« - 5e""’+"+« - 2e“’+’'"« 

_2g3+-+e 3g—+« + g—+1 + g-+2 

+e’^+“’ + 3e’“+^-« + 3e“’+^+« - e-^+^ 

0<?<lA<a; 

_2g3+-^ - 2e-’'+T + + 2e-^+'‘‘^ 

_2e’^+i^ + 2,e^+^+^ + 3e^-^ - 5e^+^-^ 

- 6""’+^ - + e“’“+2 + 2e"’+^ 


2g3+-^ - 2e-’“+T + 2e“’+^C + 26"’“+^^ 

-2e"+A + 3e-"+« + 3e"+^+« - e®+" 

_g-r. + l _ gX + 2 g-^ + 2 2e“’ + i 
_3g-+i+e _ g-+i+e 


Incidentally, this example shows also that the representation of Green’s 
operators in terms of Green’s functions—despite its long tradition in 
engineering and physics—is not always the most useful and economi¬ 
cal way of representing the Green’s operator. For many purposes it is 
better to take the Green’s operator just as an element of the operator 
ring J^4d, J] or J^[d, J^]. 
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6. Conclusion and Future Work 


While we have focused our attention to semi-inhomogeneous bound¬ 
ary problems (those with an inhomogeneous differential equation and 
homogeneous boundary conditions), one may also consider the oppo¬ 
site case of semi-homogeneous boundary problems—this is especially 
important in the case of LPDEs. The corresponding semi-homogeneous 
Green’s operator maps the prescribed boundary values to the solu¬ 
tion [IB]. In the case of LODEs, one often restricts attention to well- 
posed two-point boundary value problems (in the sense of Definition [T|). 
Writing the two evaluations as [aj and [6J and their action on u as u{a) 
and m(&), one may consider the extended evaluation matrix 


[a, h\ iu) 


/ Ui{a) 



••• 

ui{b) 

Un{h) 




which is similar to ([7|) except that it is rectangular since we consider 
more boundary functionals than we could possibly impose for one reg¬ 
ular boundary problem. If we do prescribe all 2n boundary derivatives, 
they must satisfy n relations given by the kernel of the map —)■ 

R"',X HA X ■ \_a,h\{u). For the simple example in Section [H the ex¬ 
tended evaluation matrix for the fundamental system ui = 1 ,U 2 = x 
is 


L0,ij(u) 


/l o\ 
0 1 
1 1 
VO V 


whose kernel has basis (—1,1,1, 0), (0, —1, 0,1). Written in terms of the 
boundary functionals, they encode the two relations m(1) —m(0) = m'(0 ) 
and m'(0) = m'( 1). The analogous case for LPDEs gives rise to the 
interesting notion of universal boundary problem |25j . 

There is another, more fundamental, way of extending the results 
in this paper: Currently our method for extracting Green’s functions 
works for arbitrary Stieltjes boundary problems (T, B), but only in the 
standard integro-differential algebra X = C'°°(R). Of course, it requires 
a Green’s operator G G [9, J] = X\d^ and hence a fundamental 
system for T. 
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It would be interesting to extend the concept of Green’s function 
and the corresponding extraction method to arbitrary ordinary integro- 
differential algebras 9, J). For the functional part it is 

clear how to achieve this since one sees from the structure of Green’s 
operators that necessarily g G The ring has the structure 

of a partial integro-differential algebra with derivations and integrals 

dxU ® 9 ) = {df)®g and dy{f ® g) = f ® {dg), 

rU®9) = {U)®9 and {f ® g) = f ® {\g). 

This structure will be useful for studying various properties of Green’s 
function, in particular their symmetry: For well-posed two-point bound¬ 
ary value problems (T, B) it is known [22l §3.3] that the Green’s func¬ 
tion g{x,^) is symmetric whenever {T,B) is self-adjoint. Otherwise 
one may associate to (T, B) an adjoint boundary value problem whose 
Green’s function is then g{^,x). It would be useful to know how these 
results generalizes to arbitrary Stieltjes boundary problems. 

Having an abstract integro-differential algebras the other 

problem is Green’s function will in general have a distributional part g{x,f) 
that does not £t into jF®jF. For accommodating distributions into the 
setting of integro-differential algebras, it is probably necessary to con¬ 
struct a integro-differential module generated over jF®jFhya. suitable 
notion of abstract Dirac distributions. (It is well known that distribu¬ 
tions do not enjoy a convenient ring structure, hence it seems to be 
more reasonable to go for a module. This is also the path followed in 
the algebraic analysis of T>-modules; see [SI §6.1] for example.) 
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